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FRACTURE  BEHAVIOR  OF  ICE 
IN  CHARPY  IMPACT  TESTING 


Kazuhiko  Itagaki  and  Richard  L.  Sabourin 


INTRODUCTION 

Unlike  most  structural  materials  used  in  industry, 
little  is  known  about  the  mechanical  properties  of  ice 
except  for  some  creep  and  stress-strain  data.  Several 
experiments  are  in  progress  at  CRREL  to  supply  addi¬ 
tional  data  and  to  relate  the  general  deformation  and 
fracture  behavior  of  ice  to  other  materials.  This  re* 
search  is  important  because  ice  has  the  following  pe¬ 
culiar  features:  1 )  it  is  transparent,  2)  large  single 
crystals  are  available,  3)  slip  occurs  in  most  cases  only 
on  the  basal  plane,  4)  crystals  with  relatively  few  de¬ 
fects  (104  to  10s  etch  pits/cm3)  are  available,  5)  the 
crystals  are  held  together  by  hydrogen  bonds,  6)  the 
crystals  have  a  rather  open  structure,  and  7)  a  close 
relationship  between  electrical  and  mechanical  proper¬ 
ties  has  been  observed.  Studies  of  ice  mechanics  util¬ 
izing  these  features  can  supply  useful  information  for 
the  mechanics  of  other  materials. 

Although  studies  of  mechanical  properties  of  ice 
such  as  creep,  stress-strain,  fracture  strength  and  hard¬ 
ness  measurement  have  been  made  (reviewed  by  Glen 
1975),  no  data  on  impact  testing  have  yet  been  pub¬ 
lished.  Charpy  impact  testing  is  a  standard  test  method 
for  structural  materials.  This  test  measures  energy  ab¬ 
sorbed  by  simple  beam  specimen  *  during  the  impact 
by  a  fast  single  overload  of  stress.  A  pendulum  is  used  to 
supply  energy  required  for  the  fast  stress  and  fracture 
process,  and  the  difference  between  initial  and  final 
height  of  the  pendulum  indicates  the  absorbed  energy. 
The  absorbed  energy  values  determined  are  quantita¬ 
tive  comparisons  on  a  fixed  geometry  specimen  but 
cannot  be  converted  or  generalized  into  the  energy 
values  of  other  sizes  of  specimens  or  conditions.  Gen¬ 
eral  discussions  on  Charpy  impact  testing  can  be  found 
in  ASTM  (1977). 


EXPERIMENTAL 
Sample  preparation 

Several  types  of  artificially  pown  ice  and  natural 
single  crystal  ice  from  the  Mendenhall  Glacier,  Alaska, 
were  used  for  Charpy  testing.  The  types  of  artificial 
ice  included  pure  ice,  colloidal  alumina  dispersed  ice, 
colloidal  silica  dispersed  ice,  NH4  F  doped  ice,  HF  doped 
ice,  commercial  ice  (as  received  from  a  commercial  man¬ 
ufacturer)  and  snow-ice.  All  artificial  ice  except  com¬ 
mercial  ice  and  snow-ice  were  grown  in  Lucite  tubes. 
Although  seed  crystals  were  used  to  control  the  growth 
orientation  in  these  ice  specimens,  spontaneous  nuclea- 
tion  sometimes  prevented  the  desired  result.  Most  spe¬ 
cimens  used  in  the  testing  contained  two  or  three  grains. 
No  attempt  was  made  to  orient  the  crystals  for  testing 
except  in  the  case  of  the  single  crystal  Mendenhall  Gla¬ 
cier  ice. 

The  apparatus  shown  schematically  in  Figure  1  was 
used  to  fill  the  Lucite  tube  and  to  grow  ice  crystals. 
Distilled  and  deionized  water  was  boiled  for  at  least 
one  hour  in  a  flask  for  degassing.  The  released  gas  was 
allowed  to  escape  to  the  atmosphere  through  a  gas  trap. 
Degassed  water  was  then  cooled  by  a  glass  heat  ex¬ 
changer  and  drawn  into  the  Lucite  tube  by  vacuum 
after  first  going  through  ion  exchange  resin.  A  single 
crystal  of  seed  ice  was  placed  on  the  bottom  of  the 
tube  together  with  buffer  ice.  A  small  ring-shaped 
heater  was  placed  at  the  ice/water  interface  outside 
the  Lucite  to  keep  the  shape  of  the  interface  convex 
towards  the  ice  to  prevent  sporadic  nucleation  at  the 
sail.  The  tube  was  then  lowered  into  the  cold  box  at 
the  rate  of  2.81  mm/h.  The  air  in  the  cold  box  was 
circulated  by  a  blower  while  the  air  flow  around  the 
tube  was  guided  by  ducts  to  make  a  more  uniform 
temperature  distribution.  Up  to  eight  cylinders  of  ice 
would  be  produced  simultaneously  by  this  apparatus. 
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Schematic  diagram  of  crystal-growing  apparatus.  Figure  2.  Warmed  die  Ice-sl 


The  proper  amount  of  highly  concentrated  solution 
of  eloping  materials  (HF  or  NH4  F)  was  introduced  into 
the  tube  already  filled  with  degassed,  distilled  and  de¬ 
ionized  water  to  grown  doped  ice.  Additional  evacua¬ 
tion  was  required  to  ensure  a  gas-free  solution. 

The  solutions  containing  the  desired  concentrations 
of  colloidal  silica  or  alumina  suspension  (1%  for  alum¬ 
ina  and  4%  and  1%  for  silica  by  weight)  were  boiled  for 
at  Jeast  one  hour  to  remove  the  gas.  A  system  similar 
to  that  used  for  the  pure  and  doped  ice  was  used  for 
cooling  the  solution  and  filling  the  Lucite  tube  except 
that  the  ion  exchange  column  was  not  used.  The 
growth  rate  was  the  same  for  all  three  types  of  ice. 

The  Charpy  test  specimens  were  shaped  according 
to  ASTM  standard  dimensions  (1  x  1  x  5.5  cm)  except 
for  the  notch.  Rough  cutting  of  the  notch  was  made 
by  a  band  saw,  and  then  a  warmed  die  was  used  to 
shape  a  square  cross  section  as  shown  in  Figure  2.  Etch 
pit  observation  on  this  surface  indicated  that  the  de¬ 
fect  introduced  by  band  saw  cutting  was  almost  com¬ 
pletely  removed  by  the  warm  die  and  very  few  disloca¬ 
tions  were  produced  by  this  preparation.  Uniformity 
of  specimen  size  was  another  advantage  of  this  method. 
Almost  no  attention  was  required  to  form  exact  cross- 
sectional  ice  after  proper  adjustment.  More  than  1  mm 
of  the  highly  disturbed  layer  could  be  eliminated  by 
the  die  drawing. 

Testing  procedure 

The  specimens  shaped  by  the  die  drawing  method 
were  enclosed  in  a  case  and  kept  at  a  desired  tempera¬ 
ture  for  a  prolonged  time  to  attain  thermal  equilibri¬ 
um.  The  specimens  were  then  tested  in  a  -10°C  cold 
room.  Temperature  variations  were  measured  during 
specimen  handling  and  during  actual  testing  by  a  ther¬ 
mocouple  inserted  in  a  dummy  ice  sample  having  the 
same  dimensions  as  the  test  samples.  These  variations 
were  found  to  be  within  3°,  1°,  4°,  and  6°C  respectively 
for  -2°,  -1 0°,  -32°,  and  -1 90°C  testings. 

The  specimens  were  placed  between  the  anvils  of 
the  Charpy  impact  testing  machine  (plastic  testing 
type  modified  to  accept  standard  ASTM  metal  speci¬ 
men).  The  leading  edges  of  the  specimens  were  marked 
by  a  pen  to  indicate  the  direction  of  impact  and  also  to 
make  it  possible  to  reassemble  the  shards  in  the  origin¬ 
al  order,  ft  was  recognized  after  about  200  testings 
that  a  certain  mode  of  fracture  was  related  to  extreme¬ 
ly  high  Charpy  values.  Individual  pieces  of  shards 
were  placed  in  the  original  order  and  photographed  be¬ 
fore  further  examination  for  later  measurements. 

One  of  the  central  pieces  from  each  tested  sample 
was  placed  on  a  glass  plate,  thermally  etched  and 
photom  icrographed  to  study  the  dislocation  structure 
in  the  specimen.  Concentrations  of  the  dispersed  or 


doped  materials  were  measured  using  the  remaining 
pieces  of  shards. 

Concentrations  of  dispersed  materials  in  individual 
specimens  were  measured  by  weighing  the  individual 
specimens  and  evaporated  residuals.  Only  NH4  F  and 
HF  were  used  for  the  doping  substances  and  measure¬ 
ments  of  the  fluoride  ion  concentration  were  conduct¬ 
ed.  A  fluoride  ion  electrode  (Orion  Fluoride  Ion  Ac¬ 
tivity  Electrode  94-09-00)  was  used  together  with  a 
conventional  saturated  KCI  reference  electrode.  The 
electromotive  force  between  the  electrodes  was  meas¬ 
ured  by  using  a  Beckman  model  76  pH  meter  after  a 
calibration  curve  was  made  using  known  concentration 
soiutions. 


RESULTS 
General  features 

The  histograms  representing  the  distribution  of 
Charpy  values*  of  various  types  of  ice  at  the  temperatures 
measured  are  shown  in  Figure  3,  and  the  means  and 
standard  deviations  are  in  Table  1 .  About  20  specimens 
were  used  for  each  type  of  ice  at  each  temperature.  The 
histograms  generally  show  a  normal  distribution,  though 
some  extremely  high  values  appear  in  some  cases. 

Dashed  vertical  lines  at  0.0395  ]  (0.35  In.-lbf)  serve  as 
references  for  comparison.  The  modes  of  fracture  can 
be  classified  into  three  types  depending  primarily  upon 
the  temperature  and  the  composition  of  the  specimen. 

1 .  Normal  mode.  This  mode  can  be  defined  as  a 
fracture  initiated  at  the  point  of  impact  which  pro¬ 
duces  two  shards  with  a  rather  rectangular  cross  section 
or  one  triangular  cross-sectional  shard  and  two  trape¬ 
zoidal  cross-sectional  pieces  (Fig.  4).  Charpy  test  values 
of  this  mode  generally  show  a  normal  distribution, 
though  mean  values  may  vary  depending  on  the  com¬ 
position  and  temperature. 

2.  Multiple  fracture .  Most  of  the  lowest  tempera¬ 
ture  specimens  and  some  lower  temperature  specimens 
of  HF  doped  ice  shatter  to  more  than  four  rather  rec¬ 
tangular  cross-sectional  shards.  A  typical  example  is 
shown  in  Figure  5.  The  Charpy  values  of  this  mode  of 
fracture  are  slightly  higher  and  show  a  wider  distribu¬ 
tion  than  the  normal  fracture.  Creation  of  the  new  sur¬ 
faces  may  require  additional  energy  which  must  be  sup¬ 
plied  by  the  hammer  through  certain  mechanisms.  At 
a0.01  J  (0,1  in.-lbf)  of  energy,  however,  0.1  m*  of 
fresh  surface  is  required,  assuming  that  the  energy  is 
used  only  to  create  a  new  surface  and  the  surface  energy 

•  As  the  testing  machine  Is  graduated  by  IMbf  all  readings  and 
data  procamlng  wars  made  by  In.-lbf.  Conversion  Into  SI  units 
wet  made  by  using  the  relationship  1  In.-lbf  *  0.11298  |. 


energy,  inch  pounds 

Figure  3 .  Frequency  distribution  of  Charpy  values. 


Table  1.  Mean  and  standard  deviations  of  Charpy  test  results  in  English  and  SI  units. 


Other  values  are  corrected  by  leaving  out  exceptionally  high  values  caused  by  different  modes  of  fracture. 


Figure  4.  Typical  example  of  normal  mode  fracture. 


Figure  5.  Example  of  multiple  fracture  mode. 


Point  of  impact 


Figure  6.  Example  of  fracture  from  both  ends  mode. 


is  0.1 09  J  /m2 ,  the  most  recent  value  for  ice  at  -1 0°  C 
(Ketcham  and  Hobbs  1969).  Though  this  calculation 
depends  on  the  surface  energy,  creation  of  an  unreal¬ 
istically  large  surface  is  required  to  explain  the  slightly 
higher  Charpy  values  of  this  mode,  indicating  that  the 
creation  of  new  surfaces  is  not  the  reason  for  the  high¬ 
er  Charpy  values. 

3.  Fracture  from  both  ends.  Fracture  sometimes 
does  not  start  at  the  point  of  impact  but  rather  close 
to  the  ends-presumably  aroun('  the  anvils  as  shown  in 
Figure  6.  Only  three  shards  are  produced  for  the  typi¬ 
cal  case.  The  point  of  impact  is  slightly  depressed,  and 
a  white  translucent  portion  appears  around  the  point, 
indicating  that  this  portion  has  been  highly  disturbed. 

The  Charpy  values  of  this  mode  of  fracture  range 
from  0.056  to  0.226  J  (0.5  to  2.0  in.-lbf),  which  is  2  to 
5  times  higher  than  those  from  the  other  modes.  The 
most  interesting  feature  of  this  mode  is  that  the  frac¬ 
ture  is  not  initiated  at  the  highest  stressed  point  but  at 
the  next  highest  stressed  points  at  both  ends  simultan¬ 
eously. 

Thermal  etch  pits  were  produced  on  the  fractured 


surfaces  of  pure  and  doped  ice  by  exposing  the  surface 
to  an  unsaturated  atmosphere.  Generally  the  etch  pit 
concentration  was  higher  around  the  peripheries  of  the 
fractured  surface,  the  original  surface  of  the  specimens. 
The  concentrations  were  from  104  to  10s  etch  pits/cm2 
on  the  central  portion  of  the  specimens  (within  200  pm 
from  the  original  surface),  while  on  the  periphery  of 
the  fractured  surface  the  concentrations  were  above 
the  105/cm2  range.  Dislocation  density  of  cold-worked 
surfaces  by  micro  tom  ing  or  sanding  will  be  extremely 
high  and  the  present  thermal  etch  technique  at  -10°C 
could  not  show  any  dislocation  etch  pits.  Instead,  the 
cold-worked  surface,  once  the  trace  of  cold  work  has 
disappeared  1  to  2  hours  after  microtoming,  resembles 
to  the  naked  eye  a  surface  coated  by  an  oily  substance. 
Probably  extremely  highly  concentrated  etch  pits  inter¬ 
act  with  each  other  to  prohibit  the  growth  of  individu¬ 
al  etch  pits.  Very  high  mass  transfer  along  the  surface 
will  also  play  an  important  role  in  this  phenomenon. 

The  thickness  of  this  layer  sometimes  extends  into 
the  bulk  ice  up  to  1000  pm  as  shown  in  Figure  7.  This 
photograph  is  taken  on  a  cleaved  surface  perpendicular 


Figure  7.  Etch  pits  on  the  fractured  surface  indicating  that  the  disturbed  layer 
produced  by  sanding  extends  more  than  1  mm  below  the  surface . 


Figure  8.  Etch  pits  on  the  fractured  surface  indicating  that  the  disturbed  layer 
produced  by  warm  die  shaping  is  very  shallow. 


to  the  microtomed  surface.  Warmed  die  methods  sel¬ 
dom  produce  a  high  etch  pit  density  layer  of  more  than 
200-//m  thickness  as  shown  in  Figure  8.  Etch  pits  of 
nearly  final  density  appear  a  few  minutes  after  warm 
die  drawing,  indicating  that  the  surface  destruction  on 


the  layer  is  very  small.  The  same  technique  for  reveal 
ing  dislocation  etch  pits  could  not  be  applied  to  the 
dispersed  ice  specimens  of  higher  concentrations  be¬ 
cause  segregated  particles  disturb  the  proper  etch  pit 
formation. 
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Table  2.  Linear  temperature  dependency  of  Charpy 
values  fitted  by  computer. 


Commercial  ice  C  =  0.037»(  1  - 1 .638X 1 0‘*  t) J 

|C  =  0.335  (1-1.638X  10’*t)] 

Single  crystal  (case  B)  C  ■  0.0395(1-1. 05 XI O'*  t)| 
|C  m  0.350(1-1.05X10"’t)l 

Colloidal  alumina-dispersed  C  9  0.0432(1 -1.085X  1 0'*  t)  | 
ice* _  [C  =  0.382(1 -1.085X  10~*  t)] 

♦Disregarding  the  concentration  effect. 

C -Charpy  test  value  at  t°C  in  joules. 

English  units  in  in.-lbf  are  given  in  brackets. 


Commercial  ice 

Charpy  test  results  for  commercial  ice  show  a  nar¬ 
row  symmetrical  distribution.  No  extremely  high 
values  appear  in  any  of  the  temperature  ranges  meas¬ 
ured,  although  a  slight  increase  in  Charpy  values,  up  to 
20%,  was  observed  for  the  lower  temperature  measure¬ 
ments.  No  drastic  change  was  found  in  this  tempera¬ 
ture  range,  although  the  hexagonal-cubic  transition  at 
about  -120°C  possibly  occurs  and  a  rather  linear  Charpy 
value  vs  temperature  relationship  seems  to  exist.  The 
temperature  dependency  fitted  linear  regression  is 
shown  in  Table  2.  The  narrow  distribution  of  Charpy 
values  indicates  a  uniform  structure  in  the  commercial 
ice  samples .  Therefore ,  the  effects  of  notching  and 
cold-working  by  sanding  were  studied  using  commer¬ 
cial  ice. 

Notched  commercial  ice 

A  notch  (0.4  mm  wide  x  1 .2  mm  deep)  was  cut  in 
the  center  of  individual  specimens  by  a  hot  wire  cutter 
without  introducing  excessive  defects  around  the  notch. 
The  trace,  of  the  notch  can  be  seen  on  the  left-hand 
side  of  the  specimen  in  Figure  8.  The  notch  was  faced 
opposite  the  direction  of  impact  and  tested  at  -10°C. 

A  normal  frequency  distribution  with  only  two  excep¬ 
tional  values  was  obtained,  and  the  mean  and  mode  of 
the  Charpy  values  seemed  to  agree  with  values  for  un¬ 
notched  commercial  ice. 

Two  specimens  out  of  18  did  not  fracture  from  the 
notch  while  8  out  of  18  did  not  initiate  fracture  on 
the  bottom  of  the  notch.  Whether  fracture  was  initi¬ 
ated  from  the  bottom  of  the  notch  or  from  another 
portion  did  not  seem  to  affect  the  Charpy  values. 

Sanded  commercial  ice 

The  standard  method  of  sample  preparation  in  the 
present  study  introduces  very  few,  if  any,  dislocations 
compared  to  mechanical  methods.  Dislocations  were 
deliberately  introduced  into  the  specimens  by  sanding 


with  emery  paper  (2/0)  until  the  surface  of  the  speci¬ 
men  gave  a  mat  appearance,  and  these  specimens  were 
then  tested  at  -10°C. 

The  frequency  distribution  of  the  Charpy  values  is 
rather  skewed  towards  the  lower  side,  and  the  mean 
value  is  slightly  lower  than  that  for  the  tests  on  untreat¬ 
ed  ice.  The  effect  of  the  sanding  is,  however,  quite 
small,  even  though  the  mean  appeared  to  minimal  for 
this  series  of  experiments. 

Pure  ice 

Laboratory  grown  pure  ice  was  tested  to  see  if  any 
differences  would  appear  in  comparison  with  commer¬ 
cial  ice.  Specimens  were  prepared  with  the  longer  axis 
parallel  to  the  growth  direction.  A  relatively  wide, 
slightly  skewed  distribution  with  some  scattered  extra¬ 
ordinary  values  resulted,  indicating  that  some  differ¬ 
ence  may  exist. 

Single  crystal  Ice 

All  single  crystal  ice  used  in  the  test  grew  naturally 
in  the  Mendenhall  Glacier.  Crystals  were  oriented  by 
the  reflection  of  light  from  hoar  frost  grown  on  the 
surface.  Three  combinations  of  crystallographic  orien¬ 
tation  and  direction  of  impact  were  tested. 

Case  A:  C-axis  parallel  to  the  longer  axis  of  the 
samples  (Fig.  9a). 

Case  8:  C-axis  perpendicular  to  longer  axis  of  the 
specimen  and  direction  of  impact  (Fig.  9b). 

Case  C:  C-axis  perpendicular  to  longer  axis  of  the 
specimen  and  parallel  to  direction  of  impact  (Fig.  9c). 

Case  A  shows  a  narrow  frequency  distribution  but 
skewed  to  the  lower  side.  The  mean  value  is  lower 
than  those  of  the  other  orientations  of  single  crystals 
tested.  No  extremely  high  values  appeared  except  one 
slightly  higher  value  attributable  to  the  multiple  fracture. 

The  frequency  distributions  of  cases  B  and  C  re¬ 
semble  each  other  because  the  mean  values  are  similar 
and  some  extremely  high  values  were  found  in  each 
case.  The  extreme  values  are  the  results  of  fracture 
from  both  ends  and  were  omitted  from  the  calculation 
of  mean  values  because  they  come  from  a  different 
mode. 

Initiation  of  cracks  by  pileup  of  dislocations  would 
be  easier  in  case  A  than  in  cases  B  and  C  because  the 
glide  plane  is  perpendicular  to  the  longer  axis  of  the 
specimen  and  dislocation  may  move  to  pileup  in  the 
glide  plane.  The  slightly  low  mean  Charpy  values  of 
case  A  can  be  explained  in  this  way. 

The  temperature  dependency  of  mean  Charpy  values 
of  case  B  was  studied  at  the  temperatures  of  -2°,  -10° , 
-32°  -and  -1 90°C.  The  linear  regression  fit  temperature 
dependency  is  shown  in  Table  2. 


c.  CaseC. 

Figure  9.  Crystal  orientation  and  direction  of  impact 

Snow-ice 

Snow-ice  was  made  from  fine  powder  of  ice  im¬ 
mersed  in  water  until  the  spaces  between  particles 
were  filled  by  water  and  then  frozen.  No  attempt  to 
reduce  gas  inclusion  was  made;  thus,  abundant  bubbles 
were  included  in  the  ice.  No  preferential  effect  of  bub¬ 
bles  on  the  fracture  surface  was  found.  The  mean 
Charpy  value  was  almost  the  same  as  that  of  the  com¬ 
mercial  ice.  The  mode  of  fracture  was  normal  and  in 
most  cases  two  major  pieces  were  produced. 

Colloidal  alumina-dispersed  ice 

Two  types  of  dispersed  ice  were  produced  because 
colloidal  silica  has  a  negative  surface  electric  charge 
while  colloidal  alumina  has  a  positive  charge.  Since 


colloidal  alumina  has  a  surface  electric  charge  that  is 
opposite  that  of  the  surface  of  the  growing  ice  crystal, 
easier  trapping  of  colloidal  particles  in  the  ice  can  be 
expected.  Observations  on  the  fractured  surface  after 
Charpy  testing  indicated  that  the  particles  segregated 
along  a  network  which  resembles  a  decorated  disloca¬ 
tion  network  (Fig.  10). 

Though  colloidal  alumina  disperses  in  ice  more  evenly 
than  colloidal  silica,  some  layer-like  structure  appears  in 
laboratory-grown  specimens.  The  concentration  range 
of  ice  grown  from  an  initial  1%  solution  was  from  0.25% 
to  1 .96%,  which  is  far  narrower  than  for  colloidal  silica. 
Considerable  variation  of  concentration  within  the 
specimens  can  also  be  expected  because  of  the  layer- 
like  structure  found  in  the  individual  specimens. 

The  frequency  curves  of  Charpy  values  of  alumina- 
dispersed  ice,  disregarding  concentration  variation  in 
individual  specimens,  show  a  slightly  wider  distribution 
than  those  of  commercial  ice  at  respective  temperatures. 
Linear  temperature  dependency  obtained  by  regression 
analysis  is  shown  in  Table  2  and  regression  fits  of  con¬ 
centration  vs  Charpy  values  at  each  temperature  are 
shown  in  Table  3. 

The  mode  of  fracture  was  generally  normal  with 
some  multiple  fracture  and  fracture  from  both  ends  at 
higher  temperatures.  More  than  50%  of  the  fractures 
were  multiple  at  -190°C.  Those  tendencies  are  similar 
to  those  of  commercial  ice,  though  mean  Charpy  values 
are  slightly  higher  in  colloidal  alumina-dispersed  ice. 

Colloidal  silica-dispersed  ice 

Colloidal  silica  tends  to  segregate  on  planes  parallel 
to  the  basal  plane  in  contrast  with  colloidal  alumina 
(Fig.  11).  The  concentration  of  colloidal  particles  is 
highly  dependent  on  the  growth  condition.  Distribu¬ 
tion  of  segregated  particles  seems  to  be  affected  by 
gravity  because: 

1 .  Colloidal  silica  tends  to  segregate  in  the  center  of 
the  tube  when  the  ice  is  grown  from  the  bottom  and 
the  ice/water  interface  is  concave  towards  the  ice. 

2.  Colloidal  silica  tends  to  segregate  on  the  peri¬ 
phery  when  the  ice  surface  is  convex  and  grown  from 
the  bottom. 

3.  No  colloidal  silica  is  trapped  in  the  ice  when  the 
ice  is  grown  from  the  top. 

The  concentration  of  colloidal  silica  in  individual 
specimens  is  governed  by  the  growth  conditions  and 
is  highly  variable  (0.1  to  15%).  The  colloidal  silica-dis¬ 
persed  ice  crystals  were  grown  from  suspension  with 
concentrations  of  1%  and  4%  by  weight  to  obtain  the 
respective  concentrations.  Since  the  original  solution 
(Ludox  HS-40,  DuPont)  was  supplied  as  40%  concen¬ 
trated  solution,  the  samples  were  designated  as  Vie  and 
V# o  respectively.  Presumably  the  effect  of  growth  con¬ 
ditions  prevented  this  amount  of  control  and  resulted 
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Tabte  3.  Concentration  (JxKT*)  of  dliprned  or  doped  material  vs  Charpy  value. 


Temperature  CC) 


2 

- 10 

•32 

•190 

Alumina-dtapemd  ke 

3.68+50.62 0 

4.26+49.260 

4.23+43.72D 

4. 78+50.620 

(0.326+4.480) 

|0.377+4.36D| 

(0.374+3.S7D) 

[0.432+4 .4Sd| 

SHka-dkpersed  ke 

NH4F-  doped  ice 

HF -doped  ke 

4.04+1  USD) 
(0.358+1. 14D) 

3.66+22.14D 
(0.324+1 .96D) 

4.59+4.594 
[0.406+0.4064 ) 

4.08-2.564 

(0.361-O.227d] 

4.51+S.S6D 

(0.399+0.784) 

$.92-15.704 

(0.524-1.394) 

D —concentration  in  g/mi;  d -concentration  in  mol/llter. 
English  units  in  in.-lbf  are  given  in  brackets. 


in  almost  pure  ice  for  the  1%  suspension  with  a  few 
exceptional  concentrations  of  up  to  1  JS% .  Most  of  the 
dispersed  materials  were  lost  around  the  wall  of  the 
growth  tube  or  consumed  by  the  central  part  of  the 
cylinder. 

The  frequency  curves  of  Charpy  values  for  the  4% 
colloidal  silica-dispersed  ice,  disregarding  the  variation 
of  concentration  in  each  specimen,  rather  resemble 
those  of  the  alumina-dispersed  ice.  The  mode  of  frac¬ 
ture  and  means  at  each  temperature  are  quite  close. 

The  slope  of  concentration  vs  Charpy  values  at  respec¬ 
tive  temperatures,  however,  are  about  Vs  those  of  the 
colloidal  alumina-dispersed  ice.  Differences  in  distri¬ 
bution  of  colloid  particles  in  crystals  could  have  con¬ 
tributed  to  this  lesser  strengthening. 

NH4F  doped  ice 

NH4F  and  HF  are  the  substances  known  to  make 
solid  solutions  with  ice  and  their  electrical  properties 
have  been  extensively  studied.  Hobbs  (1974)  has  com¬ 
piled  most  of  the  published  results.  For  the  mechani¬ 
cal  properties,  however,  fewer  studies  have  been  made 
except  for  some  internal  friction  measurements  (Kur- 
oiwa  1964,  Schultz  1961 ,  Walz  and  Magun  1959)  and 
creep  studies  at  lower  temperatures  (Jones  and  Glen 
1969). 

The  frequency  curve,  disregarding  the  concentration 
effect,  shows  a  wider  distribution  than  that  of  the  com¬ 
mercial  ice  with  several  exceptionally  high  values. 

Mean  and  modal  values  are  in  the  same  range  as  colloid¬ 
al  silica-  and  alumina-dispersed  ice.  Exceptionally  high 
values  are  obtained  as  a  result  of  the  mode  of  fracture 
from  both  ends.  Very  little,  or  even  slightly  negative, 
concentration  dependency  is  computed,  while  the  tem¬ 
perature  dependency  obtained  at  the  two  temperatures 
(-10°  and  -190°C)  shows  a  relatively  low  positive  coef¬ 
ficient  (about  one-half  that  of  other  types). 


Almost  all  of  the  specimens  resulted  in  multiple 
fracture  at  -190°C  while  at  -1<PC  most  specimens 
showed  normal  fracture.  A  relatively  higher  percentage 
fractured  from  bqth  ends  compared  to  other  experiments, 
however.  This  mode  of  fracture  seems  to  appear  within 
a  relatively  narrow  concentration  range  (1 .4  to  2.6x  10‘a 
mol/liter).  Six  out  of  12  specimens  within  this  range 
broke  in  this  mode,  a  considerably  higher  proportion 
than  with  the  other  types  of  ice.  A  slight  strengthening 
effect  was  observed  in  -10°C. 

HF  doped  ice 

Charpy  tests  of  HF  doped  ice  were  made  only  at 
-10°C.  An  interesting  feature  of  this  type  of  ice  is  that 
it  tends  to  shatter  into  many  pieces  (multiple  fracture). 
This  mode  of  fracture  is  observed  at  lower  temperatures 
for  other  types  of  ice.  Relatively  low  concentrations  of 
HF  seem  to  achieve  this  mode.  Eight  of  nine  specimens 
in  this  mode  were  below  a  2x  10*3  mol/liter  concentra¬ 
tion,  while  higher  concentrations  in  individual  specimens 
lay  between  0.3  and  6x  10*3  mol/liter. 

The  shape  of  the  etch  pits  on  the  fractured  surface 
of  HF  ice  is  rather  indefinite,  indicating  higher  local 
concentrations  of  HF  around  the  etch  pits.  Higher 
density  etch  pits  (well  over  10f/cma)  are  evenly  dis¬ 
tributed  on  the  fractured  surface  of  higher  (1  x  10’* 
mol/liter)  HF  concentration  specimens.  No  concentra¬ 
tion  dependency  is  found  in  this  type  of  ice,  though 
lower  concentration*  specimens  tend  to  show  multiple 
fracture,  while  specimens  higher  than  2x  1(TS  mol/liter 
in  concentration  tended  to  show  normal  fracture. 


DISCUSSION 

Three  types  of  ice  were  tested  in  a  wide  temperature 
range  (from  -2  to  -190°C).  All  showed  increasing 
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Chirpy  values  with  lowering  temperatures  within  a 
relatively  narrow  temperature  coefficient  range  (1 .05 
to  1 .638*  10*3/°C)  as  shown  in  Table  2.  No  transition 
point  was  found  within  this  range.  It  will  be  interest¬ 
ing  to  note  that  the  temperature  coefficient  of  the 
Young's  modulus  values  obtained  by  Yamaji  and  Kur- 
oiwa  (1956)  (1 .35x  10**  /*C)  and  Tabata  (1959)  (1 .395 
x  10‘3/°C)  agree  well  with  each  other  and  lie  within 
the  temperature  coefficient  of  the  present  results.* 

In  a  Charpy  test  energy  supplied  by  the  hammer  to 
the  critical  depression  s  can  be  expressed  as 


where  f  9  KEx  is  the  force  which  is  proportional  to 
depression  xtE  is  Young's  modulus,  and  AT  is  a  factor 
to  convert  the  depression  to  the  strain  in  the  specimen. 
The  Charpy  value  C  can  be  approximated  to  be  equal 
to  this  energy;  thus 


~-h 


*  ^ KEx  dx 


KEs'tl. 


0) 


This  indicates  that,  if  the  critical  depressions  remain 
constant,  Charpy  values  will  be  proportional  to  the 
Young's  modulus  which  seems  to  agree  with  the  pres¬ 
ent  results. 

The  dispersed  colloidal  materials  increase  Charpy 
values  by  the  concentration  shown  in  Table  3,  though 
this  is  not  as  large  as  those  shown  in  the  results  by 
Nayar  (1966).  The  creep  rate  measured  by  Nayar  de¬ 
creased  by  a  factor  of  30  with  a  concentration  in¬ 
crease  from  0%  to  196  of  colloidal  silica,  while  the  con¬ 
centration  range  of  the  present  study  extended  above 
10%  for  colloidal  silica  and  yet  showed  only  a  60%  in¬ 
crease  in  the  Charpy  value.  The  effect  of  colloidal  par¬ 
ticles  dispersed  in  the  ice  seems  to  differ  for  creep  and 
fracture  induced  by  impaction,  though  they  resist  the 
deformation  in  both  cases. 

There  are  some  studies  on  the  effect  of  organic  and 
inorganic  impurities  on  the  strength  of  ice.  Pounder 
(1958)  studied  the  effect  of  organic  impurities  and  ob¬ 
tained  reductions  of  strength  down  to,  in  some  cases, 
only  Vm  that  of  pure  ice.  Bailey  and  Macklin  (1 967) 


*  The  Young's  modulus  of  sea  ice  measured  by  Tabata  (1959) 
Stows  a  higher  coefficient,  especially  at  a  higher  temperature 
range.  The  high  temperature  coefficient  (5  J8x10‘*/*C)  ob¬ 
tained  by  Boyle  and  Sproud  (1931)  may  be  attributed  to  Im¬ 
pure  specimens  (grown  In  a  river). . 


obtained  leas  drastic  effects  than  Pounder  with  similar 
experiments  which  they  attributed  to  their  technique 
of  sample  preparation.  The  effect  of  HF  doping  was 
studied  for  creep  and  stress- strain  relatkmhips  by 
Jones  (1967)  and  Jones  and  Glen  (T969)  at  relatively 
low  temperatures  (-50  to  -90*C  range).  In  these  tests, 
the  creep  rate  was  found  to  be  more  than  10  times 
faster  for  ice  with  HF  doping  of  3  ppm  than  for  pure 
ice.  The  resolved  diear  stress  required  to  keep  the 
same  strain  rate  dropped  by  a  factor  of  five  for  only 
a  3-  to  4 -ppm  concentration  of  HF. 

No  such  drastic  effect  was  obtained  in  the  present 
study  during  the  Charpy  testing  of  both  HF  and  NH4F 
doped  ice.  Although  a  slight  concentration  effect  may 
result,  the  effect  is  highly  questionable  because  of  the 
rather  widely  scattered  values.  Instead,  the  effect  of 
doping  on  the  mode  of  fracture  seems  more  pronounced. 

Three  modes  of  fracture,  as  desertoed  above,  seem 
to  be  affected  by  the  composition  of  the  specimen  and 
temperature.  Most  of  the  multiple  fractures  occur  at 
lower  temperatures  except  for  HF  doped  ice.  By  using 
the  relaxation  time  data  by  Kuroiwa  (1965)  and  VanDe- 
vender  and  Itagaki  (1973),  internal  friction  peaks  of  the 
resonance  frequency  of  the  present  specimen  size 
(10,000  Hz)  can  be  calculated  since  these  peaks  should 
appear  around  -5°Qand  higher  temperature  measure¬ 
ments  (both  -2°  and  -1 0°C)  should  lie  in  the  peak  re¬ 
gion  for  pure  ice,  reducing  the  shock  wave  propagation. 
One  can  expect  lower  internal  friction  in  the  lower  tem¬ 
perature  range,  which  may  eventually  make  a  strong  in¬ 
teraction  between  the  standing  wave  or  shock  wave  pro¬ 
duced  by  impact  and  the  rather  static  loading.  Doping 
with  HF  brings  the  internal  friction  peak  down  to  -40°C 
or  lower,  and  thus  the  effect  of  low  internal  friction 
may  appear  for  the  HF  doped  ice  at  higher  temperatures. 

Values  two  to  five  times  higher  than  the  mean  re¬ 
sulted  from  the  mode  of  fracture  from  both  ends.  Ex¬ 
tremely  high  Charpy  values  are  in  most  cases  found  to 
coincide  with  this  mode.  One  end  of  the  fractured  sur¬ 
face  generally  coincided  with  the  point  of  impact  in 
normal  fracture  or  multiple  fracture,  with  an  apparently 
highly  disturbed  translucent  portion  indicating  that  the 
high  stress  under  the  hammer  plays  an  important  role 
in  the  initiation  and  propagation  of  the  crac  ;.  In  con¬ 
trast  with  these  modes,  cracks  did  not  have  any  direct 
relationship  with  the  highly  disturbed  portion  in  the 
mode  of  fracture  from  both  ends.  Crack  initiation 
seems  to  be  prevented  by  certain  mechanisms.  NH*F 
doping  in  a  certain  concentration  range  (1 .4  to  2.6x  10‘* 
mol/Uter)  or  certain  orientations  and  impact  directions 
(case  B  in  single  crystals)  seem  to  increase  the  chance 
of  occurrence  of  this  mode.  Close  investigation  of  the 
mode  of  fracture  from  both  ends  may  lead  us  to  a  way 
of  making  materials  with  higher  impact  strength. 


12 


LITERATURE  CITED 

American  Society  for  Testing  and  Materials  (1977)  1977 
Annual  Book  of  ASTM  Standards  C23  and  D256. 

Bailey,  I.H.  and  W.C,  Macklln  (1967)  The  effect  of  impurities 
on  the  mechanical  strength  of  accreted  Ice*  Journo/  of 
Atmospherk  Selene e,  vol.  24,  p.  707-710. 

Boyle,  R.W.  and  DjO.  S proud  (1931)  Velocity  of  longitudinal 
vfcratlon  in  solid  rods  (ultrasonic  method)  with  special 
reference  to  the  elasticity  of  ice.  Conod/on  Journo!  of  Re- 
search ,  vol.  5,  p.  601-618. 

Glen,  J.W.  (1975)  The  mechanics  of  ice.  CRREL  Cold  Regions 
Science  and  Engineering  Monograph  IK2b.  AD  022797. 

Hobbs,  P.V.  (1974)  Ice  physks.  Oxford:  Oxford  University 

Press 

lones,  S.J.  (1967)  Softening  of  ice  crystals  by  dissolved  fluor¬ 
ide  km*,  Physks  Letters ,  vol,  25 A,  P.  336-367. 

)ones,  S.J.  and  J.W.  Glen  (1969)  The  effect  of  dissolved  Im¬ 
purities  on  the  mechanical  properties  of  ice  crystals.  Philo¬ 
sophical  Magazine,  vol.  19,  no.  157,  p.  13-24. 

Ketcham,  W.M.  and  P.V.  Hobbs  (1969)  An  experimental  de¬ 
termination  of  the  surface  energies  of  ice.  Phf/osophkol 
Magazine,  vol.  19,  p.  1161-1173. 

Kurolwa,  D.  (1964)  Internal  friction  In  Ice.  Contributions 
from  the  Institute  of  Low  Temperature  Science ,  Hokkaido 
University ,  Ser.  A,  no.  18,  p.  1-62. 

Kurolwa,  D.  (1965)  Internal  friction  of  H^O,  0,0  and  natural 
glacier  Ice.  CRREL  Research  Report  131.  AD  615732. 

Nayar,  H.S.  (1966)  Creep  behavior  of  pure  ice  and  k*  dis¬ 
persed  with  ultraflne  amorphous  silica  (SK)t)  particles. 

Ph.D.  Thesis,  Rensselaer  Polytechnic  Institute  (unpublhhed). 

Pounder,  E.R.  (1958)  Mechanical  strength  of  Ice  frozen  from 
an  impure  melt.  Canadian  Journal  of  Research,  vol.  36, 
p.  363-370. 

Schulz,  D.  (1961 )  Die  mechanische  Relaxation  in  Eis-HF- 
Mlschkr (sullen  verscheidener  HF-Konzentratlon.  Natur- 
wissenschaften ,  vol.  48,  no.  22,  p.  691. 

Tabau,  T.  (1959)  Studies  on  mechanical  properties  of  sea  ice: 
III.  Measurement  of  elastic  modulus  by  the  lateral  vibration 
method.  Low  Temperature  Science,  Ser.  A,  vol.  18,  p.  1 15- 
129.  (Japanese  text  with  English  resume'.) 

VanDevender,  \JP.  and  K.  Itagakl  (1973)  Internal  friction  of 
single-crystal  Ice.  CRREL  Research  Report  243.  AD 
759930. 

Walz,  E.  and  S.  Magun.(1959)  Die  mechanische  Relaxation  in 
E)s-NH4F-MischkrisUllen.  Zeitschrlft  fur  Physik,  vol.  157, 
p.  266-274. 

Yamaji,  K.  and  D.  Kurolwa  (1956)  Visco-elastic  property  of 
Ice  in  the  temperature  range  0°-100*:  I.  Low  Temperature 
Science,  Ser.  A,  vol.  15,  p.  171-183.  (Japanese  text  with 
English  resuml.) 


13 


A  facsimile  catalog  card  in  Library  of  Congress  MARC 
format  is  reproduced  below. 


Itagakl,  Kazuhiko 

Fracture  behavior  of  ice  In  Charpy  Impact  testing  /  by 
K.  Itagakl  and  R.L.  Sabourln.  Hanover,  N.H. :  U.S.  Cold 
Regions  Research  and  Engineering  Laboratory;  Springfield, 

VA. :  available  from  National  Technical  Information  Service, 
1980. 

ill,  18  p.,  lllus.;  28  cm.  (  CRREL  Report  40-13.  ) 

Prepared  for  Directorate  of  Military  Programs  -  Office, 

Chief  of  Engineers  by  Corps  of  Engineers,  U.S.  Army  Cold 
Regions  Research  and  Engineering  Laboratory,  Under  DA 
Project  4A161102AT24. 

Bibliography:  p.  13. 

1.  Fracture  (mechanics) .  2.  Ice.  3.  Impact  tests. 

I.  R.L.  Sabourln,  co-author.  II.  United  States  Army  Corps 
of  Engineers.  III.  Army  Cold  Regions  Research  and  Engineering 
Laboratory,  Hanover,  N.H.  V.  Series:  CRREL  Report  80-13* 
i>u.i  smmmcnt  nwnm  o roe  .•  tm-tet-  mti 


> 

! 


